Magnetic property and crystalline electric field effect in ThCr 2 Si2-type CeNi 2 As 2 
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Millimeter sized ThCr2Si2-type CeNi2As2 single crystal was synthesized by NaAs flux method 
and its physical properties were investigated by magnetization, transport and specific heat mea- 
surements. In contrast to the previously reported CaBe2Ge2-type CeNi2As2, the ThCr2Si2-type 
CeNi2As2 is a highly anisotropic uniaxial antiferromagnet with the transition temperature Tjv=4.8 
K. A field induced spin flop transition was seen below Tat when the applied B is parallel to the c- 
axis, the magnetic easy axis, together with a huge frustration parameter / = 9w /Tm- A pronounced 
Schottky-like anomaly in specific heat was also found around 160 K, which could be attributed to the 
crystalline electric field effect with the excitation energies being fitted to Ai =325 K and A2 =520 
K, respectively. Moreover, the in-plane resistivity anisotropy and low temperature X-ray diffractions 
suggest that this compound is a rare example exhibiting a possible structure distortion induced by 
the 4/-electron magnetic frustration. 

PACS numbers: 74.70. Dd, 75.30.Gw, 75.30.Kz, 75.10.Dg, 75.20.Hr 



INTRODUCTION 



The interest in the ThCr2Si2-type structure has been 
rekindled since the discovery of superconductivity (SC) 
in (Bai_ 2; K a ;)Fe2As2i. SC was also achieved when Ba 
is replaced by other alkaline earths like Ca and Sr or 
even the divalent rare earth Eu, either by chemical dop- 
ing or pressure effect* 2 - - — . On the other hand, the nickel 
based pnictide, e.g. BaNi2As2, was reported to show SC 
too, although the nickels are non-magnetic and the T c is 
much lower—. In these 122-compounds formulated with 
ATm 2 As2 (^4=Ca, Sr, Ba or Eu, Tm=transition met- 
als), two vertically reversed Tm As layers are sandwiched 
along c-axis, while the A atoms are embedded in between, 
following a sequence of Tm As- A-Tm As. This crystalline 
structure constitutes a platform for understanding the 
interplay between Kondo interaction and Rudcrman- 
Kittel-Kasuya-Yosida (RKKY) interaction if A is re- 
placed by magnetic rare earths. Indeed, the research for 
the 4/-electron correlation in ThCr2Si2 structured com- 
pounds has been a long story, and a famous example 
is CeCu2Si2, the first heavy fermion superconductor—. 
Therefore it is very interest to study the Ce-based 122- 
nickcl pnictides like CcNi 2 As2. 

Remarkably, CeNi2As2 crystallizes in either ThCr2Si2 
(U/mmm, No.139) or CaBe 2 Ge 2 (P4/nmm, No.129) 
structure, see Fig. 1(a). The main difference between 
them comes from the interchange of the Ni and As po- 
sitions in one-half of the NiAs layers in the CaBe2Ge2- 
type CeNi2As2, which results in the loss of the inver- 
sion symmetry-ii. Previous studies on poly-crystalline 
samples have revealed that the ThCr 2 Si2-type CeNi 2 As 2 
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shows an antiferromagnetic (AFM) transition at around 5 
K, while the CaBe2Ge2- type CeNi2 AS2 is a non-magnetic 
Kondo lattice compound^ 2 -. Single crystalline samples 
are then highly desirable in order to further distinguish 
the properties of the two structures. However, since the 
occurrence of ThCr 2 Si2-type or CaBe2Ge2-type largely 
depends on the process of heat treatment^ 3 -, and in many 
cases, a mixture of the two will be derived, the properties 
of ThCi'2Si2-type CeNi2 AS2 are still not well-understood. 



In this article, we report the measurements on 
magnetic/transport /thermodynamic properties of the 
ThCr2Si2-type CeNi2As2 based on millimeter sized sin- 
gle crystalline samples. Single crystalline sample of 
CCM2AS2 was successfully synthesized by the NaAs flux 
method. We find that CeNi2As2 is a highly anisotropic 
uniaxial antiferromagnet with the transition temperature 
1V=4.8 K. The Ce 3+ moments are likely to align along 
c-axis. A field induced meta-magnetic transition (MMT) 
was seen below Tjv- Pronounced crystalline electric field 
(CEF) effect was observed. These magnetic and thermo- 
dynamic properties can be well understood by the CEF 
calculation, showing that the j = 5/2 multiplet of Ce 3+ 
splits into three Kramers doublets with the excitation 
energies Ai =325 K and A2 =520 K. In contrast to the 
CaBe2Ge2-type CeNi2As2, Kondo effect in the ThCr2Si2- 
type CeNi2As2 is not strong, with a moderately enhanced 
Sommerfeld coefficient 70 =69 mJ/mol/K 2 and a rela- 
tively low Kondo scale T% ~ 4 K. On the other hand, 
a huge frustration parameter / = 0w /?jv is obtained, 
and a frustration-distortion picture was then proposed. 
The latter highlights the important role of Ce-4/ elec- 
trons in magnetic frustrations. Therefore, the ThCr2Si2- 
type CeNi2As2 provides a new candidate for the research 
of frustration- induced magnetic and structure transitions 
and calls for more investigations. 
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II. EXPERIMENTAL 
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FIG. 1: (Color online) (a), Crystalline structure of CeNi2As2 
of ThCr2Si2-type (left) and CaBe2Ge2-type (right), (b), a 
photo of a CeNi2As2 single crystal on millimeter grid paper, 
(c), XRD structure characterization of CeNi2As2 single crys- 
tal. Only (0 21) peaks can be observed, (d), X-ray rocking 
curve through (0 4) reflection of CeNi2As2 single crystal. 



TABLE I: Lattice parameters of CeNi2As2 derived from 
Rietveld refinement based on space group U/mmm. 
a=fe=4.0806(3)A, c=9.8843(7)A, a=/3=7=90°. 
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0.86 obtained from energy-dispersive X-ray microanalysis 
(EDX) measurement. It is also comparable with previ- 
ous literature^. The derived a(b) and c are 4.0806(3)A 
and 9.8843(7)A, respectively. The small ratio c/a =2.43, 
compared to 2.80 for BaN^As^, manifests a collapsed 
ThCi'2Si2 structure. 

Quantum Design (QD) magnetic property measure- 
ment system (MPMS-5) and physical property measure- 
ment system (PPMS-9) were used in the magnetization, 
transport and specific heat measurements. Ohmic con- 
tact was made with Epoxy silver paste (Part A+B), and 
annealed in Ar atmosphere at 573 K for 30 min. Resistiv- 
ities of both in-plane (p a b, I || o,b) and out-of-plane (p c , 
I || c) configurations were measured. Thermopower was 
measured by means of steady-state technique and a pair 
of differential type E thermocouples was used to measure 
the temperature gradient. Specific heat was measured by 
heat pulse relaxation method in PPMS-9. 



High purity single crystal of ThCr 2 Si2-type CeNi 2 As 2 
was grown by NaAs flux method. CeAs, NiAs and NaAs 
were presynthesized as mentioned elsewher o 14 ' 15 . CeAs, 
NiAs, CeOz^ and Ni were weighted in the ratio of 1:1:1:1, 
thoroughly ground in an Argon filled glove box. The 
mixture was then put into a Ta tube, together with 15 
molar times of NaAs. After scaling the Ta tube by an 
Arc-melter, the tube was then sealed into a quartz tube 
filled with 0.2 bar Argon gas. The quartz tube was heated 
up to 1493 K and kept at that temperature for 10 hours, 
followed by slowly cooling down to 873 K in 10 days. 
NaAs flux was dissolved by water in a fume hood, and 
shining single crystals of CeNi 2 As 2 with a typical size of 
3 x 3 x 0.2 mm 3 were picked out from the remaining 
dreg (see Fig. 1(b)). 

The single crystalline CeNi 2 As 2 samples were checked 
by X-ray diffraction (XRD), performed on a PANalyti- 
cal X-ray diffractometer (Empyrean Series 2) with Cu- 
K Q i radiation at room temperature. Only (0 21) peaks 
can be observed (Fig. 1(c)), confirming the ThCr 2 Si 2 - 
type crystalline structure. The full width at half max- 
imum (FWHM) of (0 4) peak in the rocking scan is 
less than 0.08 °, demonstrating the goodness of sam- 
ple quality (Fig. 1(d)). We also performed Rietveld 
refinement^ on the powder XRD data (not shown). The 
structural parameters are listed in Tab. 1. The derived 
occupation of Ni site is 0.856, close to the result of 



III. RESULTS AND DISCUSSION 

The temperature dependent magnetic susceptibility 
x(T) = M/H and inverse magnetic susceptibility l/x(T) 
along B || c and B || ab arc displayed in Fig. 2(a) and (b), 
respectively. The magnitude of \ c is almost the same as 
that of Xab at 400 K, but is 16 times lager at low temper- 
ature, indicative of increasing magnetic anisotropy. Both 
X C (T) and XabiT) obey the Curie- Weiss's law above 300 
K. We fit the temperature dependent susceptibility to the 
formula x(T) = jSe — > w ^h ®w being the so-called Weiss 
temperature. The fit on the polycrystallinc averaged sus- 
ceptibility, defined as Xavg = (Xc + 2xab)/3, leads to the 
effective moment /tx e // = 2.44/ab- This value is close to 
but slightly less than 2.54 hb, the effective moment of a 
free Ce 3+ , manifesting the trivalent Ce ion and the non- 
magnetic nature of the Ni sublattice. The high magnetic 
anisotropy is also reflected in the derived Weiss tempera- 
ture, 6*^=32.2 K and 0$=-165.9 K. An AFM transition 
is signified by Xc which shows a sharp peak around 5 K 
and extrapolates to a very small magnitude in the zero 
temperature limit. Xab also shows a peak at the same 
temperature, although the reduction of Xab after the for- 
mation of the AFM ordering is much weaker. The char- 
acteristic temperature of the AFM transition, Tjv=4.8 K, 
is then determined by the peak position in the dxT/dT 
curves as shown in Fig. 2(c). We will find that this value 
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FIG. 2: (Color online) (a) Temperature dependent magnetic susceptibility of CeNi2As2 measured at B — 0.1 T shown in B || c 
(red) and B || ab (blue) directions. Poly crystalline averaged susceptibility (black) was calculated by Xavg = (Xc + 2x a &)/3. (b) 
shows inverse magnetic susceptibility. The dashed lines are guides to eyes of the Curie- Weiss fit. (c) shows the definition of 
the AFM transition temperature T N from dyT/dT, C(T) and dp/dT. (d) and (e) exhibit the evolution of AFM peak in x( T ) 
under various magnetic field, for B |[ c and B || ab, respectively. 



is also consistent with the resistivity (dp/dT) and specific 
heat (C) measurements. Fig. 2(d) and (e) show the Xc(T) 
and Xab(T) measured under various magnetic fields. It is 
interesting to notice that under increasing field the AFM 
peak shifts to lower temperatures much more faster for 
B || c than B || ab. These observations suggest that the 
Ce 3+ moments aligns along the c-axis while within the 
a&-plane the correlation between the moments is much 
stronger. The deviation of x(T) from Curie- Weiss law 
below 300 K is a sign of CEF effect and will be discussed 
later on. 

Fig. 3 shows isothermal magnetization M(B) curves 
along B || c and B || ab directions. The most fascinating 
feature for B || c is that below T/v, M(B) shows linear 
B dependence when B <2 T, and undergoes a substan- 
tial increase at around 2.4 T before a saturation trend. 
B m =2.36 T is then defined at the magnetic field where 
the increasing rate of M(B) reaches the maximum. A 
tiny hysteresis in M(B) is observed near B m , implying a 
weak first order transition. Compared with B || c, M(B) 
for B || ab is linear and much smaller, which again pro- 
vides evidence for the strong anisotropy in the magnetic 
correlation among Ce moments. This field induced MMT 
may imply the competition between the in-plane corre- 
lation and Zeeman energy. 

We now turn to the resistivity measurement. Both 
in-plane (p a b with I || ab) and out-of-plane (p c with 
I || c) resistivities were measured, and the data are shown 



in Fig. 4(a). One should notice that the ratio p c /p a b 
at 400 K is 2.3, much smaller than that of the regular 
iron pnictide A-122 compounds^^— where the ratio is 
typically of the order of 10-100. It demonstrates more 
three-dimensional electronic property in CeNi2As2, and 
is consistent with the collapsed crystalline structure. For 
T >300 K, p a b shows weak metallic behavior, while in 
T <300 K region, p a {, increases with decreasing T, and 
a broad peak centered around 110 K is observed. Sim- 
ilar behavior is also observed in p c whereas the broad 
peak position is relatively higher. Such broad peak in 
resistivity is reminiscent of the CEF effect and is further 
confirmed by the thermopower (data are not shown here) . 
The anisotropic response to the CEF in resistivities (see 
also in Fig. 9(a)) may reflect the anisotropic hybridiza- 
tion strength of electron scattering to the CEF. Another 
prominent feature is observed below 50 K, where both 
Pab and p c increase rapidly, developing the sharp peaks 
near T^r. The -logT behavior of low temperature re- 
sistivity for 5 K< T <50 K is identified in the inset of 
Fig. 4(a), revealing that CeNi2As2 belongs to a Kondo 
system with a weak Kondo scale Tk % Tjv- Fig. 4(b) and 
(c) show the isothermal in-plane resistivity p a i, versus the 
applied field perpendicular and parallel to the crystallog- 
raphy c-axis, respectively. In the case of B || ab, p a b 
decreases slightly with B at T = 10 K > T N . While 
at T — 2 K < TV, a positive magnetoresistivity (MR, 
denned as MR=(p(B) - p(0))/p(0)) is clearly exhibited. 
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FIG. 3: (Color online) (a) Field dependence of isothermal 
magnetization of CeNi2As2 in B || c, measured at 2 K, 10 K, 
and 20 K. Inset shows a tiny hysteresis near the MMT. (b) 
shows M(B) curves for both B || c and B || ab measured at 
2 K. The polycrystalline averaged M(B) is presented in (b). 



FIG. 4: (Color online) Transport properties of CeNi2As2. (a) 
Mainframe, temperature dependence of resistivity for both 
pat (I || ab) and p c (I || c). Inset displays p a b in the logT 
plot, (b) and (c), p a b as a function of magnetic field B, for 
B || ab and B || c, respectively. 



This behavior is likely associated with the suppression 
of AFM ordering under the external field. The p a b(B) 
curves for B || c are more intriguing: First, at T =10 
K, p a b decreases much faster for B || c than that for 
B || ab, providing further evidence that the magnetic 
easy axis to be c-axis. Second, at T =2 K, p a b{B) sub- 
stantially increases to a maximum near B =2.6 T, and 
then decreases drastically. The turning point of p a b(B) 
is apparently associated with the MMT observed in mag- 
netization measurement discussed previously 

In Fig. 5, we present the specific heat divided by T as 
a function of temperature. A A-shape peak is clearly seen 
at the transition temperature TV, manifesting a second- 
order phase transition. Under the magnetic field, the 
specific heat peak is suppressed to lower temperatures 
for B < B m , indicating a fingerprint of the reduction 
of the AFM transition. When B > B m , the sharp peak 
evolves with increasing magnetic field into a broad round 
peak moving to higher temperature. This Schottky-likc 
peak in C /T under a field signifies a crossover from AFM 
ordering to paramagnetic state via Brillouin-like satura- 
tion, and is consistent with the magnetic properties mea- 
surement. We also plot the C/T vs T 2 in the inset of 
Fig. 5. The Sommerfeld coefficien 70 ~ 69 mJ/mol-K 2 



is then estimated by linearly extrapolating to the zero 
temperature. This moderately enhanced Sommerfeld co- 
efficient manifests the correlation effect contributed from 
the Ce-4/ electrons. The slope of the linear extrapola- 
tion is (3 =0.498 mJ/mol/K 4 , and this leads to the Dcbye 
temperature 9d =269 K. 

In order to identify the Ce-4/ electron contribution 
to the specific heat, we consider the quantity C mag = 
Oce~(^La as shown in Fig. 6(a), where C^a is the specific 
heat of LaNi2As2 (polycrystal). It is expected that the 
difference C mag is mainly due to the magnetic contribu- 
tion because LaNi2As2 is non-magnetic and isostructural 
to CeNi2As2. As expected, a sharp specific heat jump for 
C m ag appears at T N . The jump at T N is AC mag |t=t n ~ 
6 J/mol/K. From this value the Kondo scale could be 
roughly estimated as T% ~ 4 The magnetic entropy 
gain (S m ) was calculated by integrating C mag /T over T. 
We found that S m reaches 65% of i?ln2 at X)v, and re- 
covers this value at 19 K. This suggests that the observed 
anomaly in specific heat arises from the AFM ordering of 
Ce 3+ moment in a two-degenerated ground state. Since 
Tk is low, the reduction of entropy gain at Tjv should be 
mainly attributed to the short-range ordering or corre- 
lation of the Ce 3+ moments above Tjv, rather than the 
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Kondo effect. The short-range ordering is also manifested 
by a noticeable broad tail in C mag above TV extending 
to about 20 K (see Fig. 6(a)). A broad peak centered at 
around 160 K can also be observed on C mag (T), which 
should be ascribed to the Schottky anomaly caused by 
the thermal population of CEF levels. 
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FIG. 5: (Color online) Specific heat of CeNi2As2 as a function 
of temperature, measured at various magnetic fields B || c. 
Inset displays C/T in the plot of T 2 to estimate the Sommer- 
feld coefficient 70. 

In Ce-contained compounds, the D^h (U/mmm) point 
symmetry requires a CEF Hamiltonian written as: 



V/' 

■ /(.( 



CEF — B^Ol 



B° 4 0° 4 



(1) 



where B™ {I =2,4, m =0,4) are the CEF parameters, 
while O™ are Steven's operator a 23 i 24 . In addition, the 
Zeeman interaction and exchange interaction should also 
be taken into account, 



- yf C.I 



t_L / ax ax 1 oV qV\ 
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J\\ qzcz 



(2) 
(3) 



in which x, y and z correspond to the crystallographic a, 
b and c, where g =6/7 is the Lande factor of Ce 3+ ions, 

J^ x and Jex are the components of the nearest-neighbor 
exchange interaction with Ce 3+ moment perpendicular 
and parallel to the c-axis, respectively. Combination of 
Eq. (1), (2) and (3) allows us to get the expressions of 
inverse susceptibility 2 ^, i.e., 



Xc c 



i(i + i) 



jii 
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r+iJt^Z- M-m+Q aZ). (5 , 



with j =5/2 being the total angular momentum for Ce 3+ . 
The experimental data of 1/xc an d 1/Xa& can be well 
reproduced by this model as shown in Fig. 6(b), with 



the best fitted CEF parameters B™ as well as the CEF 
energy levels and eigenstates being listed in Tab. 2. The 
ground state has a two-fold degeneracy and takes the 
form I r 7 ) = 011 I ±|) + a 2 | =f|) with <A + a\ = 

1. The calculated exchange interactions are Je X =9.l8 K 
and J e | I .=-23.1 K, respectively, manifesting anisotropic 
magnetic couplings. The broad peak in C mag (T) can 
also be well described by Schottky anomaly formula with 
energy excitations Ai=325 K and A2=520 K, and the 
result is shown in Fig. 6(a). A schematic sketch of the 
CEF split is presented in Fig. 6(d). 




FIG. 6: (Color online) (a) Magnetic specific heat C mag , de- 
rived from Cce — Clu- Solid line represents the 2- level Schot- 
tky anomaly fit based on the CEF model, where Ai =325 K 
and A2 =520 K. (b) The inverse susceptibility (open sym- 
bols) and calculated curves (solid lines), (c) Experimental 
(symbols) and calculated (solid lines) isothermal magnetiza- 
tion at 2 K and 20 K. Data for B || c are colored with red, 
while B || ab are in blue, (d) Schematic sketch of CEF energy 
levels for Ce 3+ ion in CeNi2As2. 



Above CEF analysis also allows calculating the spa- 
tial distribution of the 4/-electron charge density 2 ^ (r 7 | 
P4/(r) I T7). In Fig. 7(a), we display the iso-surface 
plot (namely the electron cloud) of the 4/-electron for 
an Ce 3+ ion surrounded by the CEF in the ThCr 2 Si2- 
type CeNi2As2. The calculation was done at T =0.1 K. 
It is evident that the electron cloud has deformed severely 
from the spherical shape, and is highly accumulated on 
the four corners which are along the Ce-As bonds, man- 
ifesting the hybridization between Ce-4/ and As-4p or- 
bitals. Another profound feature is that the Ce-4/ elec- 
tron cloud is highly "squeezed" along the z-axis (Please 
note the dimensions of x-,y-, and z-axes). It is well known 
that the topology of charge density is closely associated 
with the magnetic anisotropy. Due to the spin-orbit cou- 
pling, the orientation of the magnetic moment is coupled 
to the orientation of the 4/ charge. In the CEF the- 
ory, this anisotropy is mainly governed by the Stevens 
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TABLE II: CEF parameters, energy levels and wave functions in CeNi2As2 at zero magnetic field. 
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factor ct j 24 ' 26 . In the case of Ce 3+ ions, the negative 
Uj =-5.7143x 10~ 2 favors the magnetic easy axis to be 
parallel along the "squeezed" direction. This result rein- 
forces the previous statement that the Ce 3+ moments are 
aligned along the c-axis. Under a magnetic field B || c, 
the electron cloud is further " squeezed" along z-axis (see 
Fig. 7(b)), and thus stabilizes the original magnetic easy 
axis. In contrast, an external field B || ab elongates the 
electron cloud along the z-direction (see Fig. 7(c)), and 
consequently the magnetic moments will be rotated to 
the a6-plane. 

However, we notice that the CEF model does not per- 
fectly reproduce the isothermal magnetization measured 
at low temperatures, see in Fig. 6(c). For example, in 
the case of B || c, one expects a saturated magnetic mo- 
ment M c =6/7x5/2=2.14 //s/Ce, while the experimental 
value is about 1.36 Ms/Cc for T = 2 K and B =9 T. 
Meanwhile, we also notice that the amplitude of Weiss 
temperatures are much larger than the AFM transition 
temperature. Especially for the in-plane Weiss temper- 
ature 9$, we obtain a huge ratio f a b = 0^/Tn=34.6. 
Such a large value of f a b reminds us of the magnetic 
frustration neglected in the previous CEF analysis. A 
schematic diagram of this magnetic frustration is shown 
in Fig. 8. Taking into account only the nearest neigh- 
bor exchange interactions for both intralayer and inter- 
layer, the magnetic coupling between Cc 3+ moments are 
denoted by J\ and J2 respectively. From the CEF cal- 
culation we have a negative intralayer coupling which is 
dominating in magnitude. It means that the Ce 3+ mo- 
ments should be antiferromagnetically ordered within the 
afr-plane, as shown in Fig. 8(a). In this situation, the 
magnetic frustration stems from the J1-J2 competition 
for the moments in the two adjacent layers as within the 
extended unit cell 2 ^. We propose that a structural dis- 
tortion from the high temperature tetragonal to low tem- 
perature orthorhombic phases may possibly take place to 
release this magnetic frustration, as shown in Fig. 8(b). 
The unit cell is then doubled and the lattice constants a 
and b are no longer equivalent. This frustration- induced 
distortion scenario is reminiscent to a similar scenario in 
the iron pnictides, where the structural distortion is pos- 



sibly due to the J\-Ji magnetic frustration caused by the 
3d-electron moments 2 ^—. What we need to emphasize 
here is that the proposed frustration-induced distortion 
in the present case is caused by the Ce-4/ electrons. With 
this consideration, the CEF Hamiltonian in the low tem- 
perature orthorhombic phase can be written as 



CEF 



= b2o° 



BlOl 



BlOl 



BlO\ 



B\0\. 



(0) 



Consequently, the ground state will change into a mixed 
state like | r) =ai | ±|) + a 2 | T§) + «3 I ±5} with 
a\ + a 2 , + a|=l, and therefore a saturated magnetic mo- 
ment smaller than but close to 2.14 /is/Cc will be ex- 
pected under a moderate magnetic field. This situation 
is similar to the case in CeNiGc3 reported by Mun et 
al 3 -^. However, it is hard to get the Bp's for this new 
phase, since the information is very limited. 

To obtain more evidences for the structural distortion, 
we seek for the possible resistivity anisotropy in the ab- 
plane by performing the Pa'-pv measurement using the 
Van der Pauw's method 3 - 2 - (Please note that here a' and 
b' stand for two perpendicular directions in the afr-planc, 
not necessarily the crystallographic a and b). The mea- 
surement was carried out on a piece of square plate- 
like single crystal, and R a > and Rf were measured via 
switching the direction of electrical current. To compare 
their temperature dependences, the normalized resistiv- 
ity p' = R/R400K is used, and the result is displayed 
in Fig. 9(a). Above 200 K, p' a ,(T) and p' b ,{T) overlap 
well, while below 200 K, a discrepancy between them 
is observed. It should be pointed out that this result 
was reproduced for many times on different batches of 
single crystals. The discrepancy becomes more evident 
with decreasing temperature, exhibiting the increasing 
anisotropy in the afr-planc. It is worth emphasizing that 
this discrepancy between p' a , and p' b , can be enlarged un- 
der a moderate magnetic field B || c and reaches a maxi- 
mum near B m before it starts to decrease with further in- 
creased B (see the inset of Fig. 9(a)). Combined with the 
MMT observed in isothermal magnetization displayed in 
Fig. 3, we argue that such field dependent p' a ,-p' b , may 
be related to the field induced spin-flop transition. As is 
elucidated in Fig. 8(c), under a moderate magnetic field 
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B = B = 40 T, B//c 5 = 40 T, B//ab 




FIG. 7: (Color online) Iso-surface plot of 4/-charge density for an Ce 3+ ion surrounded by CEF in CeNi2As2. Calculated at 
T =0.1 K, (a), B = 0, (b), B =40 T, B || c, and (c), B =40 T,B || ab. The lower diagrams are the projections to the xz-plane. 




FIG. 8: (Color online) Schematic diagram of the magnetic 
structure of Ce 3+ sublattice. The adjacent Ce 3+ layers are 
denoted by colors, red and blue, while the orientations of a 
Ce 3+ moment, up and down, are signified by "•" and "x". 
The thick lines characterize the unit cell, (a) Geometry frus- 
tration stems from the J\ — J2 competition in the tetragonal 
phase, and tiny structure distortion may take place, which 
will lead to an orthorhombic phase as shown in (b). (c) Spin- 
flop transition happens when a field B || c is applied, which 
enhances the structure distortion. The arrows display the pro- 
jection of the Ce 3+ moment in the ab-plane. (d) Geometry 
frustration reduces when Ce 3+ moments are well polarized in 
the high field limit. 



B || c, the balance between the Zeeman energy and the 
AFM intralayer coupling requires that the Ce 3+ mag- 
netic moments gradually lie down to the a6-plane. This 
will then further enhance the structural distortion as well 
as the resistivity anisotropy in the afr-planc. When the 
external field is large enough, all the magnetic moments 
tend to be polarized, and thus the frustration decays with 
the increasing field. 

The frustration-distortion scenario can be further 
tested by the low temperature X-ray diffraction 
(LTXRD) experiment on the CeNi2As2 powder samples. 
We were focused on the angular range 64 < 28 < 66 °, 
where only the (2 2 0) and (116) peaks can be observed. 
The data were collected at different temperatures down 
to 12 K, the lowest temperature of our equipment. All 
the collected LTXRD patterns are displayed in Fig. 9(b). 
We find that both (2 2 0) and (1 1 6) peaks shift to the 
right-hand-side when cooling down, suggesting a shrink- 
age of the crystalline lattice, although for T <100 K, this 
shrinkage becomes very weak. We fit all these LTXRD 
patterns to a combination of two Gaussian functions, 
through which the FWHM of (2 2 0) peak is derived as 
shown in Fig. 9(c). The initial reduction of the FWHM 
with decreasing temperature should be attributed to the 
slowing down of crystalline lattice oscillation. To our 
interest, an upturn of the FWHM is clearly seen when 
T < 100 K, and especially for T <20 K, the FWHM 
increases rapidly with decreasing temperature, although 
we are not able to see the split of (2 2 0) peak directly. 
Such steep increase of the FWHM signals a tiny struc- 
tural distortion or a precursor to that happened below 
Tjv- In order to confirm the relevance of this behav- 
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ior to the 4/-electron magnetism, we also measured the 
FWHM of the LaNi 2 As 2 compound, see in Fig. 9(c). We 
find that the FWHM of LaNi2As2 drops monotonically 
with T down to the lowest temperature, in striking con- 
trast to CeNi2As2. The significant distinction betweens 
the two cases again reinforces the frustration-distortion 
possibility driven by the 4/-electrons. 




_i i i_ 
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FIG. 9: (Color online) (a), Anisotropic in-plane resistivity 
of CeNi2As2, measured in Van der Pauw's method. Inset dis- 
plays this anisotropy under various field, (b) LTXRD patterns 
of CeNi2As2 powders. The two observed peaks are indexed as 
(2 2 0) and (116). (c) Temperature dependence of FWHM 
of (2 2 0) peak. For comparison, the result of LaNi2As2 is 
also shown. 

We should remark that the structural distortion in- 
duced by the magnetic frustration of 4/-electrons is in 
general small or even tiny. One reason is that the energy 
scale for magnetic couplings is smaller than that in the 
iron pnictides. Another reason is that the next-nearest- 
neighbor coupling J3 between the Ce-moments (within 
the same layer of the Ce-sublattice as denoted in Refill) , 
which is neglected in the previous discussions, may also 
play some role. Though being small, J3 usually com- 
petes with both Ji and J2 and may lead to incommen- 
surate magnetic fluctuations^. Since a precise magnetic 
structure in CeNi2As2 is not easy to be determined by 
static magnetization and transport measurements, more 
investigations such as neutron scattering experiment are 
required to settle this issue. 

Finally, we note that CCN12P2, the counterpart com- 
pound to CCN12AS2, behaves as a typical Kondo lat- 



tice metali^. This fact demonstrates that the chemical 
"pressure", induced by replacing As with smaller isova- 
lent P 30 i 33 , promisingly acts as an effective controlling 
parameter to tune the competition between the RKKY 
interaction and Kondo coupling 3 ^. For the rare earth 
iron pnictides, such competition is much involved and 
complicated, however, mainly due to the emergent mag- 
netic order of the iron 3d-electrons2i. Owing to the ab- 
sence of magnetism in the Ni sub-lattice, the Ce-Ni based 
compounds have shown great advantage in studying the 
Ce-4/ electron correlatio n 14 ' 31 ! 35 ' 36 . This also accounts 
to the fact that the CaBe2Ge2-type CeNi2As2 is a non- 
magnetic Kondo lattice, because the Kondo coupling is 
largely enhanced by one of the inverted NiAs layer. Com- 
pared with all these cases, the ThCr2Si2-type CeNi2As2 
has a relatively small Kondo coupling but a moderately 
strong magnetic frustration. The role played by strong 
magnetic frustration on the quantum phase transition, 
in Kondo lattice in particular, remains an interesting 
issue 3 -^. Therefore, the ThCr 2 Si 2 -type CeNi 2 As 2 may 
provide a new material for the research of quantum phase 
transitions mediated by the 4/-electron magnetic frustra- 
tion. 



IV. CONCLUSION 

To conclude, we performed a systematic investigation 
on the magnetic properties and the CEF effect in the 
ThCr 2 Si2-type CeNi 2 As 2 single crystals. We find that 
this CeNi2As2 compound is a highly anisotropic uniaxial 
antiferromagnet with T/v=4.8 K. The Kondo effect is esti- 
mated to be not strong in this system, while the magnetic 
frustration of the Ce-4/ moments plays an important 
role. Pronounced CEF effect is observed in magnetic, 
transport and thermodynamic measurements. Detailed 
calculations based on the CEF theory allows captur- 
ing the electronic and magnetic properties of CCM2AS2. 
A possible frustration-induced structural distortion due 
to the Ce-4/ electrons is suggested, which is in agree- 
ment with the in plane resistivity anisotropy and low- 
temperature XRD measurements. While this issue is 
reminiscent of the frustration-induced distortion emer- 
gent in the iron pnictidc superconductors due to the d- 
electron correlation, the origin of the structural distortion 
and its relationship with the magnetic frustration in the 
4/-electron systems still need to be clarified in the future. 
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